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A model of an electrical point contact coupled to a mechanical system (oscillator) is studied 
to simulate the dephasing effect of measurement on a quantum system. The problem is solved 
at zero temperature under conditions of strong non-equilibrium in the measurement apparatus. 
For linear coupling between the oscillator and tunneling electrons, it is found that the oscillator 
dynamics becomes damped, with the effective temperature determined by the voltage drop across 
the junction. It is demonstrated that both the quantum heating and the quantum damping of the 
oscillator manifest themselves in the current-voltage characteristic of the point contact. 



There is a dramatic difference in the observed be- 
haviors of microscopic particles and of macroscopic ob- 
jects. The everyday-scale objects obey the rules of clas- 
sical Newtonian mechanics, while microscopic particles 
command the use of quantum physics for their descrip- 
tion. The effects of quantum coherence are almost never 
observed at the macroscale. The only known excep- 
tions are realized when the macroscopic quantum state 
is particularly robust against external perturbations, as 
is the case for superconductors and quantum Hall liq- 
uids. Hence, it appears natural to assume that it is the 
coupling to the external world, or environment, that 
leads to decoherence and consequently to a transition 
from quantum to classical behavior. This process was 
explored in detail in numerous works. It has been shown 
that within a phenomenological model of environment, 
at sufficiently high temperatures, a quantum mechanical 
system becomes effectively classical The environ- 

ment provides both the decoherence and the dissipation 
needed for the quantum-classical transition. 

Another important distinction between the classical 
and quantum systems is in their response to measure- 
ment. Measurement of a classical system in principle can 
have no effect on the state of the system; on the other 
hand, in the quantum regime, the measurement itself is a 
source of decoherence that inevitably changes the state of 
the system || . The main difference between the measure- 
ment process and the environment induced dephasing is 
that measurement is an intrinsically non-equilibrium pro- 
cess. In this work we demonstrate that despite the ap- 
parent differences, the measurement can also induce a 
quantum-classical transition. Recently, Gurvitz et al. Q] 
and Korotkov et al. (H| have shown that electrical mea- 
surement leads to dephasing of a finite state system that 
is being measured. The systems that they have studied 
did not have however a classical analogue. Here, we ex- 
tend their approach to the problem of the measurement 
of a mechanical system. Specifically, we consider a quan- 
tum oscillator coupled to an electrical point contact in 
the tunneling regime ||. We solve this combined prob- 
lem in the non-equilibrium limit of large voltage across 
the point contact. From the general solution, we sepa- 
rately extract the dynamics of the oscillator and of the 



current through the contact. For the oscillator, our main 
findings are that (1) coupling to the tunneling electrons 
leads to dissipation with damping coefficient independent 
of the applied voltage, (2) current shot noise generates 
fluctuations and decoherence, and (3) behavior of the 
mechanical system becomes effectively classical with the 
temperature equal to one half of the voltage drop across 
point contact. For the point contact, we determine the 
stationary non-linear current-voltage (I-V) characteris- 
tic, which implicitly measures the effects of the quan- 
tum heating and quantum dissipation in the oscillator. 
Thereby we explicitly demonstrate how measurement can 
induce a quantum-classical transition. 
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FIG. 1. Model setup. An oscillator with a confining poten- 
tial U{x) (shown as spring) is coupled to an electrical point 
contact. The motion of the oscillator modulates the electron 
tunneling between the reservoirs. The carriers in the reser- 
voirs are driven out of equilibrium by the applied bias voltage, 
eV = [iij fin > 0. 

We consider a system described by the Hamiltonian 
H = eiajai + t r a\a r + Hq(x) 

l r 

+ ^2Cl(x) (aja r + alaij . (1) 
l,r 

Here a\ (al) creates an electron at the energy level t\ (e r ) 
in the left (right) reservoir. &o(x) = — (l/2m)V| + U(x) 
is the unperturbed Hamiltonian of an oscillator of mass 
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m in classical confining potential U(x), where x is the os- 
cillator coordinate. In what follows we set both electron 
charge e and Planck constant h to unity unless stated oth- 
erwise. Generalization to multi-dimensional case when x 
is a vector is straightforward. The last term in Eq. (|l|) 
describes tunneling of electrons between reservoirs, mod- 
ulated by the oscillator position x. For simplicity we 
assume that the transition matrix element Cl (x) is inde- 
pendent of single electron states I and r in the reservoirs. 
The latter dependence can be included without qualita- 
tive modification of our principal results. Moreover, Cl(x) 
can be a general hermitian operator of the particle's co- 
ordinates, i.e. it can be a function of momenta. Finally, 
there is an electrical bias, V, applied across the tunnel 
junction, so that the chemical potentials in the reservoirs 
are related as pl — pn = eV > 0. 

To determine the time evolution of the system de- 
scribed by the above Hamiltonian, we use the many- 
body Shrodinger equation approach developed in Refs. 
for solution of highly-nonequilibrium quantum transport 
problems. The wavef unction of the full system is 

= | b (x, t)+J2 h lr{x, t)alai 

Lr 

+ W'0M)4<4' a * a *' + ■•• ( 2 ) 

1,1' r,r' 

where the initial state |0) corresponds to fully occupied 
Fermi seas in the reservoirs, with no electrons above 
the respective chemical potentials. This wavefunction 
is a superposition of all possible electron-hole combina- 
tions that can be generated by the Hamiltonian H; note 
that H conserves the total number of electrons in the 
reservoirs. The first term in \& is the time-dependent 
part of the wavefunction that corresponds to unchanged 
occupancy of Fermi seas; the second term describes a 
state in which a hole is created in the left reservoir and 
an extra electron occupies the right reservoir, etc. In 
our representation, the amplitudes b(x, t) explicitly de- 
pend on the coordinate x of the oscillator. Substitut- 
ing wavefunction (Q) into the Scrodinger equation with 
Hamiltonian ([!]), = H\ty), and performing imagi- 

nary Laplace transform on *, b(ui) — J °° dte luJt b(t), we 
obtain an infinite set of equations for the amplitudes b: 



u - H Q (x)j b (x, lS) - Q(x)t>i r (x, u) = i%(x) , (3a) 
(,r 

(uj - Ho(x) +ei - e r J bi r (x,uj) - il(x)bo(x,uj) 

— }^D,(x)bwrr'(x,uj) = 0, (3b) 



l'r' 



Here, 4>o{ x ) is the initial state of the particle. In order to 
proceed, we solve the n th equation for the n th amplitude 
b(u>, x) in terms of the (n — l) th and (n + \) th amplitudes 
and substitute it into the (n — l) th equation. Then, as it 
was shown in Ref. M, the contribution of the (n + l) th 
amplitude to the (n — l) th equation is negligible. Not 



only it is of higher order in the coupling constant SI (x), 
but it also vanishes in the limit of large bias. As an exam- 
ple, let us solve Eq. ( |3rj| ) for bi r (uj,x) and substitute the 
resulting expression into Eq. (|3a|). The sum in Eq. ( |3a| ) 
becomes proportional to 6q, 



E 



dx — — o {uj,x) , 4 

uj - h n + ei - e r 



where 4> n (x) and E n are eigenfunctions and eigenvalues 
of Hq(x), HQ(p n (x) = E n <f> n {x). Replacing summations 
over I and r by integrals over the corresponding densities 
of states, pl and pn (here assumed constant), J2i r ( tJ — 
En + ei-er)- 1 = -i7rp L p R (V + uj - E n )e{V + uj -E n ). 
The real part of the sum only weakly renormalizes the 
self energy of the oscillator M. The step- function 
can be dropped provided the typical energy spacing for 
Hq is much smaller than V and the particle's dynam- 
ics involves only relatively low- lying energy levels, such 
that condition V > E n is satisfied. The latter assump- 
tion is consistent with the final results of our calculation 
and thus Q(V + uj — E n ) is omitted. Then, using the 
completeness relation for <f> n (x), the above sum yields 
~iirp L p B {Cl 2 {x){V + uj) - Cl(x)H (x)Cl(x)]bo(u),x). Ap- 
plying this procedure to all equations in the chain (^|) and 
neglecting terms of order 0(J1 4 ) and higher, we obtain a 
set of simplified equations, 



uj - H (x)) b (x) 



In ^o,^ (z) - VCl 2 (x)Xb (x) = A(x) , (5a) 



uj — Hq{x) + ei — e r ) bir(x) — (l(x)bo(x) 



-i|{n[flb,n] {x)-vn 2 (x)Ybi r (x) = o, (5b) 
... , 

where square brackets denote commutator, 
(l[H ,Q](x) = Q(x)[H (x) x Q(x)], i] = 2i:p L p R , and 
A(x) = iipo(x)(l — it:plPr^ 2 {x)). Eqs. (||) allow us to 
derive the quantum rate equations for our system. We 
introduce the density matrix of the system defined as 

a°(x,x',i) = bo(x,t)bo(x',t), 
o^foa?',*) = ^2bi r (x,t)b* r (x',t), 

l,r 

Each of these objects is the density matrix of the os- 
cillator conditioned by the number of electrons (the su- 
perscript in cr's) that have arrived in the right reservoir. 
Here we outline the derivation of the rate equation for 
a 1 . After replacing x and uj in Eq. ( |5a| ) by x' and uj' and 
multiplying this equation by bi r (uj,x), we subtract the 



result from Eqs. (5b) multiplied by b* r (u>' , x'). Summing 



over I and r and applying the inverse Laplace transform 
to the resulting equation, we obtain 

d , 



dt 



H , a 1 



n 



?{ 

(x') - vCi 2 



(x) - FSl^a;)^ 1 
21m (Cl(x)a ' 1 ^ 
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where we introduced the off-diagonal element 
a°' l {x,x',t) = Ei,r & o( x >*)&f r 0c' ',*)• The latter can be 
expressed in terms of the diagonal element cr° . Following 
the above described procedure, we find 



,0,1 



i^{vCl(x') 



H ,n 



(x')}a° 



Substituting this expression into the equation for a 1 we 
obtain the rate equation that couples a 1 and cr° density 
matrices. The result can be generalized to the density 
matrix a n for any number n of electrons that tunneled 
across the point contact, 



a = — i 



H ,a r 



H ,n 



(x) + n 



H ,n 
h ,q 



- V (Cl 2 (x) + A 2 (a;')) }<J n - |{^(a;) 
+ Q(a/) H ,Cl (x) - 2Vh{x)h{x')^a n - 1 



(x>) 
(x') 
(0) 



From the equations of motion Eqs. (g) for the full den- 
sity matrix, one can obtain a closed form equation for the 
density matrix of the oscillator alone, a(x,x',t). Sum- 
ming Eqs. (|^) over n we find that the equation for a can 
be written in an invariant form, 



H , a 



n 



,{A,a} 



Vr] 
~2~ 



a 



(7) 



where A = [Hq, f2] and curly brackets denote anticommu- 
tator. Similarly, one can determine the time dependence 
of the average current in terms of the oscillator density 
matrix. The current is defined as (I(t)) = (N(t)), where 
(N(t)) — J2 n n I dxa n (x,x,t) is the expectation value 
of the number of electrons that have tunnelled into the 
right contact by time t. Using Eqs. (^) for & n , we find 
the current 



(/(t)> = VvjTv (0V(t)) - |Tr ([fi,A] a{t)) 



(8) 



in terms of the oscillator density matrix. Eqs. (Q) and (||) 
constitute the principal results of this work. They de- 
scribe the evolution of a system modified by its inter- 
action with non-equilibrium environment, as well as in- 
fluence of the oscillator on the current between reser- 
voirs. Remarkably, for linear coupling between the two 
subsystems, Q.(x) = f2o + Cx, we recover the Caldeira 
and Leggett equation (CL) || for the density matrix of 
the oscillator. Indeed, for the linear coupling, the sec- 
ond term in the RHS of Eq. (0) produces the dissipative 
term (r)C 2 /2m)(x — x')(d x > — d x )a, while the last term 
becomes (ijVC 2 /2)(x — x') 2 a. It is responsible for fluc- 
tuation/decohcrcncc induced by the tunnel current. 

Both the CL equation and our Eq. (Q), describe high 
temperature dynamics of a mechanical system interact- 
ing with a heat bath. There is a substantial difference, 
however. Unlike the CL work where the heat reser- 
voir is in equilibrium at high temperature T, our cal- 
culation is at T = 0. The effective temperature that 
arises in Eq. (J?]) is not the temperature of the reser- 
voirs, but rather is a result of the non- equilibrium fluc- 
tuations that arise in the course of the evolution of 



the full system described by Hamiltonian (|I|). In this 
"classical" limit, the oscillator dynamics specified by the 
master equation (Q), is given by the Langevin equation, 
mx+"fmx+\7 x U(x) = f(t), where f(t) is white noise sat- 
isfying (/(*)/(*')) = 2 ml T cS 5{t-t'). Comparing Eq. g) 
for Cl = £!o + Cx with Caldeira and Leggett equation [0, 
we deduce the damping coefficient 

7 = r]hC 2 /2m 

and the effective temperature 

T cff = eV/2. 

The fluctuations that give rise to the effective tempera- 
ture can be traced back to the shot noise of the tunnel 
current. 

In a real physical system, one can always expect pres- 
ence of additional sources of dissipation. If the current 
across the junction is the dominant source of the fluc- 
tuations for the oscillator, the effective temperature will 
remain proportional to the voltage across the junction, 
T* = T c ffj/(-f + -f r ), where 7' is the extra damping in the 
system. In this regime, the measurement of the reduced 
effective temperature combined with a separate determi- 
nation of the total damping coefficient, can provide the 
value of the non-equilibrium measurement-induced dissi- 
pation 7. 

Due to the self-consistent coupling between the oscil- 
lator and the tunnel current, the measurement-induced 
effective oscillator temperature can be directly extracted 
from the non-linear part of the I-V characteristic. The 
current dynamics is given by Eq. (||). The first term 
on the RHS of Eq. (||) is simply the Landauer formula 
with the transmission coefficient nQ 2 (x) modulated by 
the position of the oscillator. This is what one would 
expect from the Hamiltonian (|l|) assuming elastic elec- 
tron tunneling. The second term, however, is non-trivial 
and accounts for the possible inelasticity of electron tun- 
neling due to the coupling to the oscillator. For linear 
coupling fl(x) = f2o + Cx, this term yields (—7). Hence, 
the expression for the time-dependent current is 



(/(*)) 



-r l V(nl + C 2 (x 2 {t)))-e 1 . 



(9) 



For a linear oscillator with frequency loq, such that 
U(x) = (mwo/2)a; 2 , the stationary I-V characteristic can 
be determined by solving Eq. (0) for the stationary state. 
In the absence of any additional sources of dissipation, we 
obtain that in the stationary regime (x 2 ) = T c s / (mtu 2 .) . 
This generates a quadratic in V term in the I-V charac- 
teristic, which is a signature of the "quantum heating" of 
the oscillator. The constant negative offset term (—67) is 
more subtle and has to be properly interpreted. Indeed, 
for low bias this term could dominate, which would seem 
to imply a current flowing in the direction opposite to the 
applied voltage. This is however only an artifact which 
signifies the break down of our approach at low biases. 
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To trace the origin of this term we calculate the tunnel 
current in the low voltage limit at zero temperature from 
the linear response theory. Assuming that the oscillator 
is in the ground state, we obtain 

COl.r. = jV&oV + ^-(x 2 ) (eU - hwoMeV - hcu ), 

where (x 2 )o is the average square of the zero point mo- 
tion. The O-function indicates that at zero temperature, 
for small bias neither the oscillator is excited by the cur- 
rent, nor the current is affected by the oscillator's zero 
point motion. For a bias exceeding the oscillator excita- 
tion energy, an additional channel in the tunneling opens 
up. This new channel however requires the excitation 
of the oscillator, which introduces inelasticity and effec- 
tively reduces the applied voltage by Tiujq. Since the nega- 
tive offset term in the linear response equation is identical 
to the (— ej) term that appears in Eq. (^), we conclude 
that both can be attributed to the dissipative effect of 
the oscillator on the tunnel current. 

Given the correspondence between the large-voltage 
expression and the linear response result, we expect that 
the expression for current, Eq. (^|), remains approxi- 
mately valid for arbitrary voltage down to eV — Hloq. 
Below this threshold voltage, oscillator and point contact 
become decoupled and the traditional expression for the 
tunnel current in a point contact applies, assuming that 
the oscillator is initially in a ground state. In the presence 
of additional dissipation, 7', the square displacement ex- 
pectation value should be computed at the reduced effec- 
tive temperature, (x 2 ) = eVj /[2(j + j')mujQ}. With this 
modification included, the current-voltage equation (|^) 
could be used to fit experimental data in order to extract 
the measurement-induced dissipation coefficient 7. 

Based on the above discussion, it would seem impos- 
sible to measure zero point fluctuations of the oscillator. 
At low voltages, eV < Hujq, while the oscillator is in the 
ground state, it is completely decoupled from current. 
On the other hand, at high voltages the oscillator is no 
longer in the ground state. We will now demonstrate 
that it is indeed possible to measure the effect of zero 
point fluctuations using time- dependent measurement of 
the I-V characteristic. Assuming that at t — oscillator 
was in its ground state and solving the time-dependent 
density matrix equation (pi) for a linear oscillator, we find 
that the time-dependent oscillator dispersion is 

(x 2 (t)) = {x 2 ) Q(t) + {xl l (t)), 

where (x 2 )o is the oscillator ground state dispersion and 
Q(t) = e~ 7 *(sin 2 (cJo£ + <j>o) + sin 2 o;oi), with cos^o = 
7 /2w - The term (x 2 cl (t)} arises due to the "quantum 
heating." It coincides with the solution of the classi- 
cal Langevin equation with the damping coefficient rj 
and temperature T e g for initial conditions x(0) = and 
i(0) = 0, (x 2 cl (t)) = {jeV/mcj 2 .) f* dre'^ sin 2 w r for 
ujq 3> 7. In the stationary state for t = 00, we obtain 
results discussed above. For small times, t <C 7 _1 , the 



( x czW) term can be neglected and thus from Eq. (||) we 
obtain 

(I(t « 7- 1 )) = jvV(n 2 + C 2 (x 2 ) ) - e 7 . (10) 

Hence we find that the zero-point motion can be ex- 
tracted as the change in the linear part of the I-V char- 
acteristic from short to long time dynamics. 

A setup similar to the one considered here, was recently 
proposed in Ref. j^] , where the authors argued that single 
electron transistor (SET) can be used as a displacement 
detectors for a charged mechanical cantilever. It was sug- 
gested that such setup could provide the resolution suffi- 
cient to detect the cantilever zero point motion. We ex- 
pect that despite the differences in the proposed systems, 
our framework and results will be relevant for the SET 
configuration as well. Another system that may prove to 
be suitable for observation of the effects discussed here is 
a molecule trapped near a point contact 

To summarize, we have proposed a model for the cou- 
pled mechanical system and electrical point contact to 
simulate a measurement-induced quantum-classical tran- 
sition. We have solved the problem under the conditions 
of a strong non-equilibrium in the measurement appara- 
tus. Our main findings are that for linear coupling be- 
tween the mechanical oscillator and the tunneling elec- 
trons, the dynamics of the oscillator becomes damped, 
with the effective temperature determined by the voltage 
drop across the junction. We have also demonstrated how 
both the quantum heating and the quantum damping of 
the oscillator manifest themselves in the current-voltage 
characteristic of the point contact. 
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